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ABSTRACT 

The formation history of the Small Magellanic cloud (SMC) is unraveled based on the results of our 
new chemical evolution models constructed for the SMC, highlighting the observed anomaly in the 
age-metallicity relation for star clusters in the SMC. We first propose that evidence of a major merger 
is imprinted in the age-metallicity relation as a dip in [Fe/H]. Our models predict that the major 
merger with a mass ratio of 1:1 to 1:4 occurred at ~7.5 Gyr ago, with a good reproduction of the 
abundance distribution function of field stars in the SMC. Furthermore, our models predict a relatively 
large scatter in [Mg/Fe] for — 1.4 < [Fc/H] — 1.1 as a reflection of a looping feature resulting from the 
temporally inverse progress of chemical enrichment, which can be tested against future observational 
results. Given that the observed velocity dispersion (~ 30 km s" 1 ) of the SMC is much smaller than 
that (~ 160 km s _1 ) of the Galactic halo, our finding strongly implies that the predicted merger event 
happened in a small group environment that was far from the Galaxy and contained a number of 
small gas-rich dwarfs comparable to the SMC. This theoretical view is extensively discussed in the 
framework that considers a connection with the formation history of the Large Magellanic cloud. 
Subject headings: Magellanic Clouds — galaxies: formation — galaxies: star clusters — stars: abun- 
dances 



1. INTRODUCTION 

Merging plays a key role in the hierarchical galaxy 
formation in the cold dark matter (CD M) Universe 
(jWhite fc Reeslll978t IWhite fe Frenklll991l ). Numerical 
simulations based on the CDM scenario dem onstrate how 
the galaxies build up hie rarchically (e.g., lAbadi et al.l 
120031: ISpringel et al J 120051) . In addition to the theoret- 
ical ground, the observed mass function of galaxies as 
well as that for cluster of galaxies have evolved with red- 
shifts, suggestive of the merging process over the cos- 
mic time, though some discrepancy between the predic- 
tion and the observation exists (e.g..lPozzetti et alJl2007t 
iMarchesini et al J 120081: iBahcall fc Bodell2003D . Turther- 
more, some fraction of galaxies at all redshifts show an 
ongoing merger. For high-z galaxies, we witness an asym- 
metric structure composed of many b lobs, indicative of 
a major merger (jConselice et al.ll2003D . In the near field, 
recent minor merger can be seen as a tidal stream in 
the halo such as the Sagittarius stream in the Galaxy 



(Ibataetal 



(Ibata et al 



20011) . 



1994) or the Giant Southern Stream in M31 



Recent studies have suggested that Magellanic clouds 
(MCs) are the pro mising candidate s to ex hibit a sign of 
past merger. First. [Bekki fc Chi ba (2008|) have claimed 
that kinematical difference be tween H I gas and older 
stars in the Small MCJSMC) (jStanimirovic et al.ll200l 
lHarris fc Zaritskvl 120061) is an end result of a major 
merger occurred in the early stage of the SMC forma- 
tion. Secondly, the Large MC (LMC) is found to possess 
an extended stellar halo exhibiting a broad metallicity 
spread, which can be linked with the nature of hierarch i- 
cal formation through mergers (|Majewski et all l2009h . 
Finally, the LMC and S MC might have belong ed to a 
larger Magellanic group p'Onghia fe Lakdl200l . which 
ensures a gas infall onto the MCs owing to the smaller ve- 



locity dispersion inside the MC group. Note that such an 
infall is hard to realize once the MCs become the Galaxy 
satellites like the present time. 

Chemical features of MCs have been extensively inves- 
tigated through the obserations for young and old stars 
belongi ng to the fi eld as well as star clusters (see the re- 
view bv lHillll200i ). though there are less studies for the 
SMC where its past information can be obtained only 
through the star clusters. The merit of clusters as a tool 
to trace the chemical evolution is that their ages esti- 
mated from the color-magnitude diagram are far more 
precise than those of field stars. Accordingly, how the 
metallicity evolved with time, i.e., the age-metallicity re- 
lation of star clusters, deserves the first-ranked informa- 
tion for the study on the history of MCs. The past rav- 
aging event such as a major merger, if occurred in galax- 
ies, will imprint some signature in their age-metallicity 
relation. 

In the general scheme on chemical evolution of galaxies, 
the metallicity gradually increases with time. If we spot a 
break in this relation for a galaxy, what does this feature 
tell us about the formation history of the galaxy? In this 
letter, we report that there exists a clear dip in [Fe/H] 
in the age-metallicity relation among the star clusters 
in the SMC, and conclude that it is a relic of a major 
merger occurred at the early stage of SMC formation. 
Its implication for the history of SMC in connection with 
the LMC and the Galaxy is also discussed. 

2. EVIDENCE AND MODEL 

2.1. A dip in [Fe/H] 

Filled circles in Figure 1 denote the correlation of 
[Fe/H] with age for individual star clusters. Here we 
have compiled the latest data from various references, 
and the result for 386 field stars that arc divided into six 
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Fig. 1. — Comparison of age-metallicity relation between the 
models and the observations. Three curves represent the results 
of the model in which the merger with the mass ratio of one to 
four is included (solid curve) , the model with an equal mass merger 
(dashed curve), and the model with no merger event (dotted curve), 
respectively. Filled circles denote the observed data of star clusters, 
compiled from various papers (Da Costa & Hatzidimitriou 1998; 
Piatti et al. 2001, 2005, 2007a,b; Kayser et al. 2007; Glatt et al. 
2008a, b), whereas open circle s show the mean met allicity in six age 
bins for 386 field SMC stars llCarrera et al.ll2008aT) . 
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Fig. 2. — Age-metallicity relation for old (t > 5 Gyr) clus- 
ters by three different groups. Black circles represent the data 
of Glatt's group extracted from Fig. 1 for the clusters with their 
ages older than 5 Gyr. Red an d bl ue circles are t aken from 
IDa Costa fe Hatzidimitrioul [IT998T) and IMiehell et all 1T99SI) , re- 
spectively. The overall trend implied by the age-binned averages 
(crosses; t < 7 Gyr, 7 < t < 10 Gyr, t > 10 Gyr) is shown by the 
arrows. Here we neglect three data situated around (t, [Fe/H])=(6, 
-1.1) (see the text). 



age bins is also shown (open circles) . A fairly good agree- 
ment between clusters and field stars is seen. As an over- 
all feature, we see a gradually increasing trend of [Fe/H] 
with time. At the same time, a gap in level between 
the evolutionary path for t >8 Gyr and that for t <7 
Gyr is conspicuous. In other words, the Fe abundance 
once decreases ~8 Gyr ago and then restarts to increase 
around t—1 Gyr till the present. Here we want to stress 
that this dip is not an artifact resulting from the system- 
atic difference between the multiple studies. Except for 
one c luster (at t—3.5 Gyr, IDa Costa fe Hatzidimitrioul 
fl99cl . ages and metallicities of other clusters are the re- 
sults of two groups by Glatt's et al. and Piatti et al., 
each of which occupies the separate age-range in the 
plot, i.e., for old (t > 6 Gyr) clusters, and young and 
intermediate-age clusters, respectively. Thus, the ob- 
served dip is established by one study. Moreover, the 
previous results on old clusters are also suggestive of the 
presence of dip. Figure 2 shows the data of Glatt's group 
extracted from Fig. 1 for the clusters with their ages older 
than 5 Gyr (black circles), t ogeth er with the results by 
Da Costa fe Hatzidimitrioul (|1998l ) (red circles) and by 
Migh elTet al.l (|1998f ) (blue circles). We see the dip con- 
vincingly for one while marginally for the other. Mean- 
while, this figure makes the inhomogcneity of metallicity 
values around t ~ 6-6.5 Gyr by Glatt's group stand out. 
The lower value of [Fe/H] for one cluster seems compat- 
ible with other two studies. Finally, it should be of note 
that the dip feature is supported by the data of field stars 
as well as by an increasing rate of [Fc/H] for young and 
intermediate-age clusters as shown in §3. 

Only a mechanism to make such a dip in [Fe/H] is 
considered to be an external dilution of the interstellar 
matter (ISM) by an accretion of very low-metallicity gas. 
The degree of dilution is not small: the decrease in [Fc/H] 
is estimated to be 0.3 dex at least. That is equivalent to 
an accretion of pristine gas with an equal mass to the 



existing ISM at that time, likely an order of 10 8 M Q , as 
the least mass scale. Therefore, the observed decrease in 
[Fe/H] in the SMC has a strong implication for the past 
occurrence of a major (not minor) merger. 

In the following, we construct the model for chemical 
evolution of the SMC, based on the major merger sce- 
nario. 

2.2. Model 

Our prime aim is to reproduce the age-metallicity re- 
lation of the SMC, by incorporating the major merger 
phenomenon into the model. Let ip(t) be the star for- 
mation rate and A(t), A' if) be the gas infall rate, the 
accretion rate by a major merger, respectively, then the 
gas fraction f g it) and the abundance of heavy-element i 
Ziit) in the gas change with time according to 

^k = -am+A{t)+A'it) (1) 



= -aZ.it^it) + Z A>i {t)A{t) + yu,ii>{t) 

+Vu,i I i>{t-t la )g(t la )dt 1& + Z AI>i {t)A'(t) , (2) 
Jo 

where a is the mass fraction locked up in dead stellar 
remnant and long-lived stars, yi is the heavy-element 
yield from an SNII or SN la, and Za,% (Za'j) is the 
abundance of heavy element contained in the infalling 
gas (the accreting gas driven by a merger). 

The star formation rate tp{t) is assumed to be propor- 
tional to the gas fraction with a constant rate coefficient 
of 0.03 Gyr" 1 . For the infall rate, we adopt the formula 
that is proportional to exp(— t/T 1B ) with a timescale of in- 
fall Ti n . In this formula, T- m — 0.3 Gyr is adopted, and the 
metallicity Za,% of an infall is assumed to be [Fe/H] =-2 
with a SN-II like enhanced [a/Fe] ratio. For reference, 
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Fig. 3. — Observed and predicted abundance distribution func- 
tions of SMC stars against the Fe abundance. The observed dis- 
tribution represented by crosses with error bars is obta ined from 
the selection of data taken from Carrera ct al. (2008aj). See the 
text for detail. The model distribution denoted by filled circles is 
the one calculated under the assumption of the merging with the 
mass ratio of one to four. The other one (open circles) is with the 
one-to-one merger model. These distributions are convolved using 
the Gaussian with a dispersion of 0.1 dex in [Fe/H], and are nor- 
malized to coincide with the total number of the sample stars of 
observation. 



the values for the corresponding rate and timescale in the 
solar neighborhoo d are predicted t o be 0.4 Gyr -1 and 5 
Gyr, respectively (|Tsuiimotoll2007t) . 

The effect of a merger on chemical evolution is devised 
to be an additional term A'(t) which has the same form 
as A(t). Here, the onset of A'{t) is set at t = 5.5 Gyr 
with the assumption that the age of SMC is 13 Gyr old. 
Thus, A'(t) = for t < 5.5 Gyr, A'(t) oc exp(-(i - 
5.5)/rin) for t > 5.5 Gyr. The timescale of accretion is 
assumed to be r; n = 1 Gyr. This somewhat long timescale 
is adopted to avoid the sudden decrease in [Fe/H] by a 
merger. The total gas accretion is normalized so that 
J Q (A(t) + A'(t))dt = 1. The mass ratio between two 
functions is the crucial factor to control the impact of 
the merging. We calculate three cases of A : A' = 1 : 
0, 1 : 1, 1 : 4. In all cases, the metallicity Za' is set to be 
[Fe/H] =-2 so as to produce the dip in [Fe/H] by a merger 
event. 

For the initial mass function, we assume a power-law 
mass spectrum with a slope of -1.35, which is combined 
with the nucleosynthe s is yie lds of SNe la and II taken 
from iTsujimoto et al.l (]1995l ) to deduce ynj and J/i a ,i- 
We apply an instantaneous recycling approximation to 
SN II progenitors, leading to the ejection rate of heavy 
elements, which is proportional to ip(t). In contrast, the 
effect of time delay yields the ejection rate from SN la's, 
which is proportional to tp(t — t\ a ). It is assumed that ti a 
spans over some range according to a distribution func- 
tion g{t\a). Here we adopt that the fraction of the stars 
that eventually produce SNe la for 3 — 8M& is 0.05 with 
a box -shaped g(tj a ) for 0.5 < t/ a < 3 Gyr (|Yoshii et al.l 
Il99fih . 

3. RESULTS 
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Fig. 4.— Predicted correlations of [Mg/Fe] with [Fe/H] for SMC 
stars, together with the observed data at present (filled circle: Venn 
1999, open circle: Rolleston et al. 2003). Solid and dashed curves 
are the same as for Fig. 1. 



First, we show the calculated age- metallicity relations 
for three cases (dotted curve: no merger, dashed curve: 
one-to-one merger, solid curve: one-to-four merger), to- 
gether with the observations. Here ages represent the 
lookback time. The merger models can reproduce not 
only the observed dip around 7 Gyr ago but also the in- 
creasing rate of [Fe/H] before and after the dip. There 
are several clusters deviating from the evolutionary paths 
predicted by the merger models. These are classified into 
two groups; one is well fitted by the model with no merger 
(two clusters) and the other is composed of ones dropping 
below the curves due to their lower metallicities than the 
predictions at each age (the assemblage of four clusters 
at t ~ 2 Gyr and one at t ~ 6.5 Gyr). The presence 
of such clusters might imply that there existed regions 
where a gas accretion driven by a major merger did not 
take place, and that additional accretion events, possi- 
bly by minor mergers, locally diluted the ISM in the 
SMC, along with another possibility that low- metallicity 
clusters came into the SMC from other galaxies through 
major/minor mergers. 

The successful merger models should be verified by 
checking whether they are compatible with other chemi- 
cal features of SMC. In general, stellar abundance distri- 
bution function (ADF) gives a stringent constraint on 
the models. For the SMC, from over 350 field stars 
(jCarrera et al.l 12008 aT ). the observed ADF can be con- 
structed. The distribution denoted by crosses in Figure 
2 shows the derived ADF form their data, in which we ex- 
clude stars in three outer fields where the metallicity on 
average is distinctly lower than other fields. In the end, 
293 stars are entered. The ADF thus obtained is com- 
pared with the results by the one-to-four merging model 
(filled circles) and the one-to-one model (open circles). 
A fairly good agreement with the observation, especially 
for the one-to-four model, can be seen. 

Finally, we show the resultant correlation of [Mg/Fe] 
with [Fe/H] for two merging models. After drawing a 
loop which is a result of the increase-decrease-increase 
chain of [Fe/H] induced by a merger, both curves reach 
the present-day observed value. For the pr esent value , 
chemical compositio ns for A-type supergian t (|Vennll 19991 ) 
and B-type dwarf (Rolle ston et al.l I2003D are shown. 
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Our merging models predict a relatively large scatter in 
[Mg/Fe] for the range of -1.4 < [Fe/H] < -1.1 as a reflec- 
tion of a looping feature. 

4. DISCUSSION AND CONCLUSIONS 

Motivated by the observed dip in [Fe/H] seen in the 
age-metallicity relation for the star clusters in the SMC, 
we construct the chemical evolution model in which the 
gas infall driven by a major merger is incorporated. Ac- 
cordingly, we claim that the SMC experienced a major 
merger in the manner that two progenitor galaxies of the 
SMC merged with the mass ratio of 1:1 to 1:4 around 7.5 
Gyr ago. Different speed of chemical enrichment between 
two galaxies results in creating the gap/discontinuity 
of age - [Fe/H] relation as observed. lHarris fc Zaritskvl 
(2004) found that the SMC formed relatively few stars 
between 8.4 and 3 Gyr ago and thus did not claim 
the presence of a starburst population around 7.5 Gyr 
ago. This together with our finding implies that a ma- 
jor merger occurred in the SMC did not trigger a major 
starburst for some physical reasons but proceed with a 
moderate star formation. We do not doubt that the de- 
tailed processes involved in SMC evolution, including the 
connection between the potential merger events and the 
star formation history, will be unraveled by the future 
databases of detailed elemental abundances of old stel- 
lar population in the SMC, which will be brought by the 
new- generation of surveys such as HERMES. 

Here we have proposed a galaxy-galaxy merger in the 
framework that the SMC merged with a gas-rich dwarf 
galaxy including a stellar component, since such process 
would change the dynamical properties of the SMC, the 
signature of which is suggested to be imprinted in the 
kinematical difference between H I gas and old stars 
(jBekki fc Chibal l2008h . If we put aside the kinemati- 
cal discussion, an accretion event of a massive HI gas 
cloud with a mass of 10 8 M Q order, which is not asso- 
ciated with a galaxy embedded in a dark matter halo, 
also explains the observed [Fe/H] dip in the SMC. How- 
ever, such a massive gas accretion event is almost out- 
side the range of possibility. That's firstly because such 
massive isolated gas clouds are extremely rare and usu- 
ally associate d with tidal debris as a result of interacting 
galaxies (e.g.,[Koribalski ct al. 20041). Nowhere in the Lo- 
cal Group has possessed a potentiality to provide such a 
huge amount of gas with a short timescale of ~ lGyr. 
On the other hand, if the SMC was already within the 
Galactic halo about 7.5 Gyr ago, a gas accretion onto the 
SMC rather than onto the Galaxy would be severely sup- 
pressed owing to the large relative-velocity between the 
SMC and the cloud in the Galactic halo. In the end, it is 
concluded that the observed dip is attributable to the gas 
infall associated with a galaxy-galaxy merger occurred to 
the SMC. 

The difficulty of a gas accretion onto the SMC in- 
side the Galaxy discussed above also holds for a merger 
case. Gi ven that the observed ve locity dispersion (~ 30 
km s"\ lHarris fe Zaritskvl 120061 ) of the SMC is much 
smaller than that (~ 16 km s _1 estimated from y2V c , 
iBinnev fe Trem ainc 1983) of the Galactic halo, the major 
merger is highly unlikely to happen within the Galactic 
halo. Accordingly, we suggest the following two possi- 
ble schemes on the SMC's merger event. One is that the 
merger event could have happened in a small group where 



(i) the velocity dispersion was comparable to that of the 
SMC and (ii) there was a number of gas-rich dwarfs. 
Recent cosmological simulations have shown that the 
Galaxy grows th rough accretion of small groups (e.g., 
iLi fc Helmil I2008T ). which is compatible with the view 
that the small group including the SMC might have ex- 
isted before and accreted into the Galaxy potential to 
lose other member galaxies in the Galactic halo. The 
other is that the SMC had been a binary system of two 
gas-rich dwarfs until ^7.5 Gyr ago and then merged with 
each other to form a single dwarf with an extended HI 
gas disk. 

By incorporating the hypothesis that the MCs have 
a common diffuse dark halo with the mass larg er than 
~_2x 10 10 M Q into the orbital models for the MCs. lBekkil 
(2008) succeeded in reproducing the present 3D velocities 
of MCs, which are c onsistent with the latest proper mo- 
tion measurements (jKallivavalil et al. Il2006f) . With the 
same model, he found that the strong tidal interactions 
among the LMC, the SMC, and the Galaxy for the last 
2 Gyr led to the formation of the Magellanic Stream. It 
is, however, unclear how the possible common halo was 
formed in the histories of the LMC and the SMC. As the 
possible mechanism to form the common halo, it can be 
regarded as the remnant structure evolved from a group 
where both the LMC and the SMC had belonged until 
recently. This scenario leads to a plausible picture that 
such a group provided the SMC with the major merger 
event as proposed in the present study. We note that 
the recent dynamical coupling of the LMC and the SMC 
around 3-4 Gyr ago is another promising mechanism for 
the formation of the common halo in terms of the origin 
of "th e age gap" in the globular cluster system of the 
LMC (jBekki et al. Il200l . 

Interestingly, we see a similar dip in [Fe/H] at ~8- 9 
Gyr ago for field stars in the LMC (|Carrera et al. 2008b), 
though in contrast to the continuous age-distribution 
of field stars, the LMC star clusters show the lack of 
the intermediate-age (4 < t < 9 Gyr ) clusters (e.g., 
lOlszewski et~aLlll991l: iDirsch et al.ll2000t ). Whether this 
observed dip in [Fe/H] for the LMC has some connec- 
tion with that for the SMC at ^7.5 Gyr ago is the is- 
sue to be solved with a high priority. If the LMC's 
dip is due to a merger event occurred at the slightly 
older age than the SMC, the event will imprint some 
observable records in stellar structure and kinematics of 
the present-day LMC. The potential candidates are con- 
sidered to be the ob served extended stellar halo (e.g., 
iMaiewski et al . 2009) a nd the pr esence of the thick disk 
(e.g.. Ivan der Marel et al."l l2002| ). In the end, we can 
bear the following picture in mind. A small group, so- 
called the Magellanic group, contained both the LMC 
and the SMC at least until the completion of a merg- 
ing process about 6.5 Gyr ago, and then it had been 
partly or completely disrupted by the Galactic tidal field 
during its accretion onto the Galaxy. It results in the 
present situation that the other group member galaxies 
exist as satellite gal axies with orbits similar to the MCs 
in the Local Group (|D'Onghia fc Lakell2008t) . Validation 
of this scheme together with the issue of which satellite 
galaxies originally belonged to the Magellanic group will 
be discussed in our future papers. 
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